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Abstract 

By using an ansatz for the eigenf unction, we have obtained the exact an- 
alytical solutions of the radial Schrodinger equation for the pseudoharmonic 
and Kratzer potentials in two dimensions. The energy levels of all the bound 
states are easily calculated from this eigenfunction ansatz. The normalized 
wavefunctions are also obtained. 
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1 Introduction 

It is well known that the study of exactly solvable problems has attracted much 
attention since the early study of non-relativistic quantum mechanics. The exact 
solutions of the Schrodinger equation for a hydrogen atom (Coulombic) and for a 
harmonic oscillator in three dimensions as well as in an arbitrary number of spatial 
D-dimensions represent two typical examples in quantum mechanics [1-4]. The 
pseudoharmonic and Mie-type potentials [5,6] are also two exactly solvable po- 
tentials other than the Coulombic and anharmonic oscillator. In general, there are 
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a few main methods to study the exact solutions of quantum mechanical systems. 
The first traditional method is to solve the second order differential Schrodinger 
equation by transforming it into some well known ordinary differential equations, 
whose solutions are the special functions like the confluent hypergeometric func- 
tions, associated Laguerre polynomials, the Whittaker functions and others [2]. 
The second method is related with the SUSYQM method, further closely with 
the factorization method [7]. The exact quantization rule also shown its power 
in calculating the energy levels of all the bound states for some exactly solvable 
quantum systems such as, the Morse, the Rosen-Morse, the Kratzer, the harmonic 
oscillator, the first and second Poschl-Teller and the pseudoharmonic oscillator 
potentials [8-10]. The Nikiforov-Uvarov method [11] which is introduced for the 
solution of Schrodinger equation, to its energy levels, by transforming it into hy- 
pergeometric type second order differential equations. The method is based on 
the determination of the solution of wavefunction in terms of special orthogonal 
functions for any general second-order differential equation [12-15]. Also, the 
perturbative method introduced for calculating the energy levels of some bound 
states of the Schrodinger equation for some approximately solvable quantum sys- 
tems [16-19]. Further, the exact solution of the £)-dimensional Schrodinger equa- 
tion with some anharmonic, pseudoharmonic and Kratzer potentials are obtained 
by applying an ansatz to the wavefunction [20,21]. 

Ozfelik and §im§ek have proposed an ansatz for eigenfunction and presented 
the exact analytical solution of radial Schrodinger equation in three-dimensions 
of all the bound states with inverse-power potentials [22] . Moreover, by applying 
this factorization ansatz to the eigenfunction, Dong has carried out an exact ana- 
lytical solutions of the Schrodinger equation in two-dimensions for some physical 
inverse-power potentials [23]. 

It is well known in quantum mechanics, a total wave function provides im- 
plicitly all relevant information about the behaviour of a physical system. Hence 
if it is exactly solvable for a given potential, the wave function can describe such 
a system completely. Until now, many efforts have been made to solve the sta- 
tionary Schrodinger equation with anharmonic potentials in three dimensions and 
two dimensions. With the same spirit, in this work, we apply the method to study 
the solutions of the Schrodinger equation in two-dimensions for some pseudo- 
harmonic and the modified Kratzer potentials due to their applications in physics 
[5,6,8,10,21,24]. 

This paper is organized as follows. In Section [2l we apply the factorization 
ansatz to the eigenfunction to obtain an exact analytic solution to the stationary 
radial Schrodinger equation in two-dimensions for the molecular pseudoharmonic 
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and modified Kratzer-Fues potentials. We also obtain the bound state eigensolu- 
tions of these molecular diatomic potentials by making a suitable ansatz to every 
wave function. The concluding remarks will be given in Section [3l 



2 Wavefunction Ansatz 



Consider the two-dimensional stationary Schrodinger equation with a potential 
V{p) depending only on the distance p from the origin [25]: 



2fi 



V' + V{p)]ij{p,v) = Eij{p,^), 



(1) 



leads to the following equation for the radial part of the wave function [25,26]: 

^(p,<^)=0. (2) 



+ -— + -^ + ^{E-V{p)) 



dp"^ p dp p^ difp- Ti 
On expressing the wave function as the product 

i){p, if) = r~^^'^Rm{p) exp{±irrnp), m = 0, 1, 2 
the resulting differential equation for Rm{p) reads 

(m^ - 1 /4) 



cP 2p 
_dp^ h 



— + -^{E- V{p)) 



Rm{p) = 0, 



(3) 



(4) 



where m and E denotes the angular momentum and the energy spectra, respec- 
tively. From the works [22,23], the radial wavefunction Rm{p) has the following 
ansatz: 



where 



/m(p) 



Rm{p) = /m(p)exp [g{p)] 



1 when m = 0, 

IlT=i(P - when m = 1, 2, 3, 



and 



g{p) = ap^ + 61np, a < 0. 
On inserting Eq.(5) into Eq.(4), one obtains 



(5) 



(6) 



(7) 
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9"{p) + ig'ip))' + 



fM + V'm{p)9'{p) 

fm{p) 



Rm{p) = 0, (8) 



where the prime denotes the derivative with respect to the variable p. Following 
the Refs [27,28] it is found from Eq. (4) that R'^{p) can be expressed as 



fl fl P^ Jm{P) 

which is the most fundamental equation for the following analysis. 



2.1 The pseudoharmonic potential 

This potential has been studied in three dimensions [5] using the polynomial so- 
lution method and in D-dimensions using an ansatz for the wave function [21]. It 
has the following form [5,21]: 

V{p)^dJ^-^) , (10) 

\Pe P J 

which can be simply rewritten in the form of isotropic harmonic oscillator plus 
inverse quadratic potential [5,10,29] as 

y(r) = V + 4 + C, AS > 0, (11) 
p 

where A = D^p-^, B = Depl and C = -2De. 

First of all, we study the ground state (m = 0) so that we take /o(p) = 1 and 
g{p) in Eq. (7) to solve Eq. (9) 



C'-E' + A'p' + ^ = 2a (1 + 26) + ia'p' + (12) 

p"^ p^ 

where E' = ^E, A' = ^A, B' = ^B and C = ^C. On equating coeffi- 
cients of pP, in Eq. (12), for all p = —2, 0, 2, the relations between the potential 
parameters and the coefficients a and b are as follows: 

B' -^^b'^ -b, A' = 4a^ -E' + C = 2a(l + 26). (13) 
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Hence, we find from the first equation in (13) that b — ^ + y/W and from the 



second equation, a 



'-, but only the negative root yields a regular wave 



function at p = 0, and of course we choose a — so that A' > 0. Finally, 

from the third equation in (13), the energy is 



£;o = c + 



/i 



1 + 



(14) 



and the corresponding eigenfunction, by virtue of Eq. (3), can be written as 



1 

^2 P 



(15) 



Secondly, for the first node (m = 1), we use = p — ci{^ and g{p) given in 
Eq. (7) to solve Eq. (9), 

C'-E' + J^p" + = 2a (3 + 26) + 4aV' + ^i^±ll. (16) 

p'^ p^ 

The relations between the potential parameters and the coefficients a, h and a^^^ 
are 

+ - = 6^ + 5^ ^' = 4^2^ -E' + C' = 2a(3 + 26), ^ = 0, (17) 
which are solved as 



6 = -^ + V5' + l,a = - 
So that the energy eigenvalue reads 



A', = C" + /^(3 + 26). 



(18) 



/2rA 



P 



£;i = c + 

and the corresponding eigenfunction can be written as (cf. (3)) 



(19) 



1 \2pA 2 

-2 P 



(20) 
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We follow the same steps for the second node (m = 2) (i.e. /2(p) — {p—a['){p— 
af'')) and g{p) is defined in Eq.(7) which imply that the following algebraic equa- 
tion holds 



R' 4- — 

C'-E' + A'p^+ " 



, =2a(5 + 2t) + 4«y+ "'+^f + '\ (21) 



which leads to equations 
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B' + ^ = (6 + l)(6 + 2), y4' = 4a^ - E' + C ^ 2a{5 + 2b) , af^ = ^ = 0. 

(22) 

Therefore, the energy eigenvalue is 

^2 = C+ w 1 1 + 

and the eigenfunction reads 



2/i5 



+ 4 , 



(23) 



1 /2_M 2 
^2 ^ 



(24) 



It is now possible to understand the general structure of the calculation with arbi- 
trary value of the angular momentum m. One has the general algebraic equation 



R' _L ^2 _ 1 



C'-E'+A'p'+ 



which yields the equations: 



= 2a (2m + 2b+ l)+4a>M 



(6 + m) (6 + m — 1) 



(25) 



B' + m'^-- = {b + m){b + m-l), A' = Aa^, -E' + C ^ 2a{2m + 2b+l). 

(26) 

Finally, the general formula for energy eigenvalue reads 



Em — C + 



'2h^A 



1 + 



2p.B 



+ m2 , m = 0,1,2, ••• 



(27) 



and the corresponding normalized radial wave function is factorized in the form 



1 rijlA 



2V r 



(28) 



where all normalization constants can be evaluated from the condition 



rdr — 1. 



(29) 



For example, the explicit calculation shows that 



N 



2J^) 



which is the normalization constant. 



(30) 



2.2 The Mie-type potentials 

This potential has been studied in the D dimensions using the polynomial solu- 
tion and the ansatz wave function method [6,21]. An example on this type of 
potentials is the standard Morse [30] or Kratzer-Fues [31,32] potential of the form 
[6,21,24,33] 

Vip) = -De - f ) , (31) 

where D^, is the dissociation energy between two atoms in a solid and is the 
equilibrium internuclear seperation. The standard Kratzer potential is modified 
by adding a term to the potential. A new type of this potential is the modified 
Kratzer-type of molecular potential [21,24] 

Vip) = D,(^^^y , (32) 

and hence it is shifted in amount of D^. The potential in Eq. (32) has been stud- 
ied in D dimensions [21] by making the wave function ansatz [20]. However, 
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this potential has also been discussed before in three dimensions [24] and in D 
dimensions [6]. This potential [33] can be simply taken as 



A B 
V{t) = - + - + C, 
P P 



(33) 



where A = -D^Ve, B = D^rl and C = [21,24]. 

First of all, for the ground state (m = 0), we take /o(p) = 1 and g{p) in Eq. 
(7) to solve Eq. (9) 



p p 



2 2ah 6(6-1) 



(34) 



P P 

where E' = ^E, A' = ^A, B' = ^B and C = '^C. On equating coefficients 
of for all p = —2,0,2, the relations between the potential parameters and the 
coefficients a and h are 



B' --^h{h- 1), A! = 2a6, - E' + C ^ o?. 



(35) 



Hence, we find from the first equation in (35) that 6 = | + \f^' and from the 



second equation in (35) is solved by a 
reads 



A' 



1+2VB' 



=. Finally, the energy eigenvalue 



2fiAyh^ 



Eq — C — 



and the normalized wave function is 



1 + 2 



2ij,B 



2' 



V'o(P> = -^0 exp 



2/x 



-^{E- C)p 



I 2fj.B 



(36) 



(37) 



Secondly, for the first node (m = 1), we use — p — cti^'' and g{p) given in 
Eq. (7) to solve Eq. (9), 



(38) 



Hence, the relations between the potential parameters and the coefficients a, h and 

(1) 
a\ ' are 
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+ - = 6^ + 6, A' = 2a(l + 6), -E' + C' ^ a^ ai^ = 0, 
and so the exact energy eigenvalue is 



Ei^C- 



1 + 2a/^ + 1 



2 ' 



(39) 



(40) 



and the corresponding wave function is 



2/i 



i^{E- C)p 



,/^+ig±i<^_ (41) 



We follow the same steps for the second node, m — 2 (i.e. /2 (p) = (p — ^ ) (p 
q;2^^)) and g{p) is defined in Eq.(7) to find 



p p^ 



and their solutions are hence found to be 
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+ ^ = (6 + 1) (& + 2) , A' = 2a{h + 2), - + C" = a^, af^ = = 0. 

(43) 

Therefore, the energy eigenvalue is 

E2 — C — 



2f,A'/Tr 



1 + 2^/^+4 



2 ' 



(44) 



and hence the corresponding wave function is 

V^2(P,<^) = iV2exp 



2p 



pV a^ +^e±*2^. 



(45) 



Thus, in general, for any arbitrary m, we find that 



„ ^, , A\ B' + m'^ - I _ ^ ^ 2a{b + m) , {b + m) {b + m - 1) 



C'-E' + — + 
P 



+ 



, (46) 
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which leads to the following algebraic equations 



1 



B' + m^-- = {b + m){b + m-l), A' ^ 2a{b + m), -E' + C ^ a\ (47) 



Finally, the energy eigenvalue is 

Em — C — 



1 + 2, /^ + m2 



2 ' 



(48) 



and the corresponding normalized wave function is 



where the normalization constant 



2/x 



P 



(49) 



N ^ 



{2^-^^{E-C)) 



(50) 



2j^ + m2 + l)! 



3 Concluding Remarks 

We have easily obtained the exact bound state solutions of the two-dimensional 
radial Schrodinger equation for two general potential forms representing the pseu- 
doharmonic [5,21] and modified Kratzer molecular [6,21] potentials by using the 
wave function ansatz method [22,23]. The presented procedure in this study is 
systematical and efficient in finding the exact energy spectra and corresponding 
wave functions of the Schrodinger equation for any desired quantum system. This 
method is simple in producing the exact bound state solution for further anhar- 
monic potentials. 
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